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We calculate the charge transport in a flux biased dissipative Cooper pair pump using the method 
of full counting statistics (PCS). This is instead of a more traditional technique of integrating a very 
small expectation value of the instantaneous current over the pumping period. We show that the 
rotating wave approximation (RWA), which fails in the traditional technique, produces accurate 
results within the PCS method. 
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I. INTRODUCTION 

Mesoscopic quantum systems evolving periodically and 
adiabatically in time provide an opport unit y to study 
geometric eflFects, e.g., the Berry phas^^l^ in a dis- 
sipative environment. In particular, normaP^ and 
superconducting charge pumping devices belong to this 
class of systems. The relation between the pumped 
charge and the Berry phase was proved for supercon- 
ductingpumps by Aunola et al.^^. In experiments with 
qubits^ and atoms^, in order to observe the Berry phase, 
superpositions of energy eigenstates must be generated. 
In contrast, in the charge pumping devices the system 
remains all the time close to the ground state. What 
is, then, being measured is the variation of the ground 
state's Berry phase with the change of the phase bias 
of the system. The first experimental determination of 
the Berry phase in a superconducting Cooper pair pump 
(CPP) was performed by Mottonen et aZ.™ 

The adiabatically controlled CPPs are susceptible to 
environmental no ise an d dissipation. These effects have 
been investigatecP^^. One common approach to de- 
scribe the influence of dissipation in such systems is to 
employ the quantum master equation technique. This 
way one obtains the time evolution of the reduced den- 
sity matrix of the system. However employing RWA in 
order to obtain a Markovian master equation in Lindblad 
fornpSl^ which grants complete posi tivity , leads to incor- 
rect results for the pumped ch&vg^^^. To obtain the 
correct pumped charge one has to go beyond the RWAES, 
i.e. include the fast rotating terms of the interaction pic- 
ture master equation. The very simple reason for this 
is the fact that the tiny instantaneous value of the cur- 
rent is determined by the tiny off-diagonal elements of 
the density matrix dropped within RWA. 

In this paper we employ an alternative formalism of 
FC^^. This is achieved by including (theoretically) an 
ideal measuring apparatus in our description of the elec- 
trical circuit, i.e., a capacitor with infinite capacitance. 
Such a capacitor has zero impedance at all frequencies, 
i.e., it does not disturb the circuit's dynamics. Yet it 
allows to introduce formally the number of Cooper pairs 



that have passed through the system and the conjugate 
counting field. It turns out, that within this approach the 
RWA can be used safely and we obtain a master equa- 
tion for the combined system of CPP and the measuring 
device in a Lindblad form. Thus, we develop a formal- 
ism which results simultaneously in the correct pumped 
charge and a positive definite reduced density matrix of 
the system. We use the results of Ref. [T3]as a benchmark 
(obtained by going beyond the RWA) and obtain a per- 
fect matching with our current results. In this context 
it should be mentioned that in superconducting phase 
biased devices the interpretation of the higher FCS cu- 
mulants is tricky-^''. 

After providing a brief introduction to the CPP system 
in Sec.|TTj in Sec. |III| we introduce a gedanken measuring 
device, which allows us to use the methods of FCS. We 
briefly demonstrate the power of this approach by repro- 
ducing the connection between the Berry phase and the 
pumped charge for a closed system^. In Sec. IV we derive 
a Markovian quantum master equation for the reduced 
density matrix of the combined system consisting of CPP 
and the measuring device. It should be emphasized that 
in the derivation the RWA was used with respect to the 
fast rotating terms in the interaction picture. In Sec. [V] 
we present and discuss the obtained results and compare 
them with previous calculations, which have been done 
via integrating the currenl!^. In Section VI we conclude. 



II. THE SYSTEM 



The sy stem we are studying is the so called Cooper 
pair sluic^iSliSUSI depicted in Fig. [l] (a). It consists of 
two SQUIDs enclosing a superconducting island. The 
island is controlled via a gate voltage Vg whereas the 
SQUIDs are manipulated by external fluxes <I>i. The sys- 
tem is well studied theoreticalljB^ and was also realized 
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FIG. 1. (a) Flux biased Cooper pair pump consisting of two 
SQUIDs enclosing a superconducting island. The SQUIDs 
are controllable by external fluxes 4>i and c&cxt. The island 
is driven by a gate voltage Vg. (b) Evolution of parameters 
Ej\, Ej2 and rig during a pumping cycle of period r = 80 ps. 
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FIG. 2. Cooper pair pump as in Fig.[T] An infinite capacitor 
C„i has been implemented as a measuring device for trans- 
ferred charge creating a new phase (pm to the system. 



experimentalljfl^. The Hamiltonian is given bjff^ 



H = Ec{n - iigf - E,ii{(t>i) cos [ + ip 



~2~ 



i?J,2(</>2)cOs( ^-^). (2.1) 



Here, h is the operator of number of excess Cooper 
pairs on the island, Ug = CgVg is the gate charge, 
Ec ~ (2e^)/(C) denotes the charging energy with C be- 
ing the capacitance of the island, Ejj{(j)i) are the tun- 
able Josephson energies of the SQUIDs, 0oxt is the total 
phase of the circuit and (p — {ipi — ip2) corresponds to 
the phase on the superconducting island. The connection 
between the fluxes and phases is given by (j>i = 27r<i>i/$o 
with $0 being the magnetic flux quantum. Since the 
phase If and the number operator n are conjugated vari- 



ables, the commutation relation 



''^ holds. 



To pump Cooper pairs through the device, the external 
parameters Eji and Ug have to be altered cyclically and 
adiabatically. An example of a pumping cycle is given 
m Fig.[l](b). In the be ginning of the cycle the SQUIDs 
are closed and there is no gate voltage applied to the 
island. First, the left SQUID is opened by increasing 
its Josephson energy to the value of Ey^^. After the 
left SQUID is opened, the gate charge Ug is increased to 
pull a Cooper pair onto the island. In the next step the 
left SQUID is shut off and the right SQUID is opened. 
Decreasing the gate charge finally pushes the Cooper pair 
into the right lead, so that altogether one Cooper Pair is 
pumped through the device. 



Real SQUIDs are not perfectly closable, therefore we 
consider a residual Josephson energy i?™™ = 0.03Ey^^, 
which is a realistic experimental value. This causes a 
supercurrent flowing through the island in addition to 
the pumped current. 



III. IMPLEMENTATION OF THE MEASURING 
DEVICE 

To calculate the transferred charge, we use the tech- 
nique of full counting statisticJ^^ (PCS). In particular, 
we consider a gedanken experiment in which we add a 
charge measuring device to the circuit. This device is 
a capacitor Cm described by a phase variable lys^ (see 
Pig. [2|. To minimize the disturbance of the circuit by 
the measuring device, we choose Cm ~> oo, such that no 
charging energy term is added to H. 

To demonstrate the functionality of the capacitor as a 
measuring device, we have a closer look at the (reduced) 
density matrix pM of the measuring device. Whenever 
N Cooper pairs reach the capacitor, its state changes 
as |</Jm) — >■ e.^'^'^'^\fm) ■ Hence the off-diagonal element 
\'~Pm){'^'m\ of the density matrix changes as 

I \/ / I N Cooper pairs transported iNr(m -^m' \\ \l I \ 

(3.1) 

If we now assume a probability distribution P/v of iV 
transferred charges, the corresponding transformation of 
the density matrix of the measuring device reads 

N 



where x{X) is the desired cumulant generating function 
(CGP) to calculate the transferred charge. Por the count- 
ing fleld we get A = ipm — f'm- analogy to Ref. [HI the 
CGP is obtained as 

Trp [{ip 

m 



X(A) 



(3.3) 



Trp [(<p™|pp(0)|O]- 

With the CGP, the cumulants of the number of Cooper 
pairs transferred onto the capacitor Cm are calculated as 

C„ = (l/z)"ann(x(A))U=o. (3.4) 

In particular, we are interested in the charge Q — 2eCi. 
After the incorporation of the capacitor Cm into the sys- 
tem, the new Hamiltonian for the system plus measuring 
device reads 
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H 



PM 



dipmHp{(prn) {(fim) {(fir. 



(3.5) 



with 



Hpi(pm) = Ec{n - Ugf - Ej^i{(j)i) cos 



The property of HpM being diagonal in ip^ 
absence of the conjugate charge is due to - 
means that Lp,n is a constant of motion. 
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(3.6) 
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and rig g [0, 1], it is possible 
to restrict the Hamiltonian to a two level Hilbert space of 
two charge stated {|0), |1)}, where |0)(|1)) correspond 
to zero (one) excess Cooper pair on the island. The two- 
level Hamiltonian reads 



where we assume an adiabatic ground state evolution, i.e. 
\^^{T,<Pm)) - e-*®°('^-)e-»QB(..™)|g(o,(p^)). Here &u = 
/g dtEg{ipmj t) is the dynamical phase of the ground state 
of the system and 9b = —§ {g{q)\'^ ^\g{q))dq is the ge- 
ometric Berry phase of the ground state^. The vector q 
belongs to the reduced parameter space spanned by two 
parameters q = (a^(<y9m), 7(Vm)) ■ From the first cumu- 
lant we find for the transferred charge 



= -2eCi = 2e (d^^QD + d^„M 



(3.13) 



We identify the first term of eq. (3.13) as the charge 
Qg = 2e9(^^^OD due to the super-current. Because the 
pumped charge solely depends on the time dependence of 
the system, the second term of eq. ( 3.13 ) Qp = 2ed^^^ 9b 
corresponds to the pumped charge. This concurs with the 
previous result^. 



IV. LINDBLAD MASTER EQUATION 



H 



2ls 



-1/2 



Ecil-2ng)a, 



[Ejl + Ej2) cos I z I fJrr 



{Eji ~ E,n) sm ( | ct„ 



(3.7) 



By diagonalizing this Hamiltonian we get the energies 
Egj^ for the corresponding ground- and excited states 



|<7(¥'m))=ae^^|0)+6|l), 
|e(^™))=6e'-'|0)-a|l). 



(3.8) 
(3.9) 



Here, the amplitudes a and h as well as the phase 7 are 
real numbers satisfying + 6^ = 1, which depend on the 
controlled parameters of the Hamiltonian {Eji,Ej2,ng) 
as well as on the "measuring" phase ipm- They are ex- 
plicitly given by 



(ng - 1/2) 



{ng-l/2f + §i- 



7 = arctan 



E 



■72 



E 
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E 
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E 



tan 



J2 



Ooxt - <y5,, 



(3.10) 



(3.11) 



where we introduced the abbreviation E12 — 
y/Ej^+ Ej2 + 2EjiEj2Cos{(t)c^t- fm)- A similar 
representation without the measuring device was used 
in Ref. US 

As a first verification of the technique of FCS, we cal- 
culate the transferred charge in the closed system, i.e. 
without coupling to an environment. We consider one 
pumping cycle of period r. It follows from eq. (3.3 1, 
that the CGF is given by 

;^(A) = e-^®°('^'")e"'®^('^"'e*®=^(^'")e'®°(^"\ (3.12) 



In this section, we derive a Markovian master equa- 
tion of the Lindblad form to describe the influence of 
dissipation on our system. The derivation is performed 
along the lines of Ref. [TH a detailed calculation including 
the RWA can be found in appendix |X] The main source 
of dissipation are fluctuations SVg of the gate voltage 
Vg = Vg + SVg, which results in a coupling proportional 
to az in the charge basis. These are, e.g., due to the dis- 
sipative parts of the controlling circuit impedance. Also 
material dielectric losses can be effectively modeled this 
way. In previous calculationJ^^ the secular approxima- 
tion (RWA), well known from the field of quantum optics, 
turned out to be inadequate (non charge conserving)!^. 

Our starting point is a Markovian master equation in 
the interaction pictur^l^ 



dt 



PPKlit) 



dsTiB[Hi{t), [Hi{t~s),ppM{t)®PB]] ■ 

(4.1) 



Here, ppuit) is the density matrix of the combined 
system of pump (P) plus measuring device (M) and 
Hj{t) — ApM{t) B{t) describes the interaction be- 
tween the system, represented by ApA,/(i)and the bath, 
represented by B{t). The explicit form of ApM in the 
Schrodinger picture is 



(4.2) 



ApM = / dlfrnEc<Jz <^ | </?,„)( V^m | • 



The density matrix pp of the bath is assumed to corre- 
spond to a thermal equilibrium state. 

To perform the RWA we need to consider eq. (4.1) in 
an explicit basis. Because of the cyclic time dependence 
of the control parameters with period r, we use the Flo- 
quet states \(f>n{t)) for this purpose. These states are 
cyclic solutions (up to a phase) of the Schrodinger equa- 
tion with the property |0„(i -|- r)) cx |0„(i)). They can 
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be obtained numerically and have the useful propertjEl^ 
that they implicitly contain su per adiabatic corrections 
to the instantaneous eigenstatea^il. For a review on Flo- 
quet theory, see Ref. HOI 



In order to compute the CGF according to eq. (3.3) 



it is necessary to evaluate the time evolution of the re- 
duced density matrix of the Cooper pair pump includ- 



ing the off-diagonal elements with respect to the mea- 
suring (counting) phases {ipra\ppM{t)W„i) = /Cp"'^" (*)■ 
In the limit Cm oo, the counting fields ipm and (p'^ 
are constants of motions. This means that the time evo- 
lution factorizes, i.e. Pp^'^'^it) depends only on matrix 
elements with the same Lpm and (p^. This significantly 
simplifies the analysis of the master equation as we only 
have to deal with 2x2 matrices. We obtain 



H'f + H' 



LS 



7(0) 



T ^mHJ„^ t! \ I T \ T n^" 

^oPp ^0 ^ 2 V " 



pT^'-l','l'' 



E 

j#fce{e,g} 



7, ^jkPp ^jk 7^ ^jkPp ^jk 



^ML],L,kPT^'- + 7(^.)pr^"i;-fc^i;- 



jk ^jk 



(4.3) 



Here, 7(w) and ^(w) arc the real- and imaginary 
parts of the half-sided Fourier (Laplace) transform 
of the bath correlation function 7(a;)/2 -|- i£,{uj) — 
{B{s)B{0)) e^^^ ds. The transition frequencies are 
given by ujjk ~ Ek{t) — Ej{t) [the explicit form is given 
in eq. (A5)], and we introduced the Lindblad operators 
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Lo{Vm,t) ^ Ec^{(j)n\(Tz\(j)n)\(l)n){<i)n\, (4.4) 

n 

Ljk{'Pm,t) ^ Ec{cj)j\(Jz\(t)k)\cl>,){(likl (4.5) 
as well as the Lamb shift 

(4.6) 

Quantities labeled with and without prime depend on ip'^ 
and Lpm respectively. The Lindblad operator Lo(ipm,t) 
describes the pure dephasing, whereas Ljk{ipm,t) cause 
relaxation and excitation. The coherent part of the mas- 
ter equation consists of a "commutator" between pp""^™ 
and Hp'' + Hls (this is not exactly a commutator since 
the counting phase in the Hamiltonian iJp " -I- Hf^s is dif- 
ferent depending on which side of the density matrix it 
is placed). 

It should be emphasized that in, contrast to Ref. |I2] 
and[T31 we did perform the secular approximation (RWA) 
concerning fast rotating terms in the interaction picture. 
The resulting master equation for cpm = f'm is therefore 
a Lindblad master equation. For ip^ ^ ^^e matrix 

Pp'"'^"'(t) has neither trace of unity nor is it hermitian, 
therefore, a Lindblad form cannot be expected for the 
off-diagonal entries. 



A typical description of the gate voltage fluctua- 
tion is provided by an ohmic bath, where the spec- 
tral density is linear in frequency. Furthermore, we as- 
sume the bath to be in thermal equilibrium 7(0;) = 

7ow ^1 — exp • For comparison with previ- 

ous worlPS the duration of a pumping cycle is assumed 
to be r = 80 ps, the charging energy of the island is 
Ec/{27rh) = 21GIIz with the ratio between the maximal 
Josephson energy of the SQUID and the charging energy 
given by Ej^^^/Ec = 0.1. The external phase is chosen 
to be (jjcxt — — 7r/2. In this parameter regime the two 
level description as well as adiabaticity are well justified. 
The energy splitting of the system along the pumping 
contour is shown in Fig. |4] (a). Temperatures of order of 
the minimal energysplitting uj^in ~ 2 GHz x 27rfi should 
be sufficient to significantly influence the system. 



We now go on to numerically solve the master equa- 
tion. In order compare with experiments, where multiple 
pumping cycles are performed, we allow the system to 
reach the quasistationary stat^^. To do so, we start 
pumping in the ground state. After a certain number 
N of pumping cycles {N depends on temperature and 
the coupling strength to the dissipative bath) the system 
approaches a quasi-stationary state in which the density 
matrix is periodic, pp{t + r) = pp(t). We then investi- 
gate the charge transferred during time t after the quasi- 
stationary state is reached. It is obtained from equations 
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(3.3) and (3.4) and is given by 



QWtot = 2e^-ln 

OA I 



Trp 




\ 




Trp 




I 


A=0 



(5.1) 



where we numerically approximate the differentiation 
with respect to A by a finite difference quotient. 

The pumped charge can be distinguished from the su- 
percurrent by the fact that the former changes sign when 
the pumping cycle is traversed in the opposite direction, 
while the latter does not. The pumped charge was calcu- 
lated as Qp = (Qtot - Qtot)/2, where Qtot is the charge 
transferred when the pumping cycle is traversed in the 
opposite direction. It is plotted as a function of time 
within one pumping cycle in Fig. [3] (a). 

To understand the behavior of the pumped charge as 
a function of the coupling strength Eq^q to the gate 
voltage fluctuations we analyze energy splitting during 
the pumping cycle (Figj4^ (a) ) and the matrix element 
\{(j)e{t)\cr^\(j)g{t))\ (Fig. |4p(b)). In a time independent 
situation the relaxation rate would be proportional to 
\{(j)e{t)\(7z\4'g{t))\'^ . Upon inspection of Fig. Wl (a) and (b) 
it is obvious that the influence of dissipation reaches a 
maximum, whenever the energy splitting is minimal. On 
the other hand, whenever the energy splitting is large, 
the relaxation rate is low and the system has not enough 
time to relax. This indicates that a stronger coupling to 
the environment should support ground state pumpinj^^l. 
We note that the pumped current is largest at times when 
the voltage parameter Ug crosses the charging energy de- 
generacy point at rig = 0.5. The reason is that at these 
times the state of the system changes at the highest rate. 
Furthermore, we see that decreasing the coupling to the 
environment leads to a smaller pumped charge. This can 
be explained by the fact that with a low coupling to the 
environment the system gets easily excited during the 
times with low energy splitting. At times with high en- 
ergy splitting, the relaxation rate becomes too small for 
the excitation to relax. After several cycles the quasi- 
stationary population of the excited state can therefore 
be quite large. As the excited state carries charge in the 
opposite direction the total pumped charge is reduced. 

A similar trend is also visible in the total charge trans- 
ferred through the island as depicted in Fig.[3](b). We 
see that a supercurrent flows through the island at all 
times. It reaches its maximum a.t t — 20, 60ps, where the 
gate voltage is Ug = 0.5. Comparing both figures in detail 
we observe that the transferred charge due to the super 
current is also influenced by the dissipative environment. 
This is apparent, as the decrease of the pumped charge is 
smaller than the decrease of the total transferred charge. 

In Fig. |5] (a) the pumped charge Qp{t) for one com- 
plete pumping cycle is plotted as a function of the cou- 
pling strength for different temperatures. As mentioned 
before, for weaker coupling strengths as well as for higher 
temperatures the pumped charge tends to decrease. As 




FIG. 3. (a) Pumped charge Qp through the island during one 
pumping cycle for fixed temperature ksT j {2nh) — 2 GHz and 
for different coupling strengths Eq'Jo — 0.2(solid), Eq^o = 
O.l(dashed) and Eqjo ~ 0.05(dotted) plotted over time. As 
expected, charge is only pumped when the gate voltage is 
changed. The weaker the coupling to the environment, the 
more the pumped charge tends to decrease, (b) The total 
charge Qtot transferred through the island. The temperature 
is kpT / {2Tih) = 2 GHz and it is plotted for different cou- 
pling strengths Eq^q = 0.2(solid), -Bc77o = O.l(dashed) and 
E'cIo ~ 0.05(dotted). The total charge also tends to decrease 
as the coupling strength is decreased. 




FIG. 4. (a) Energy splitting of the two level approxima- 
tion of the system for one pumping cycle, (b) The quantity 
\{4'e['t)WA4>g{'t))\ is plotted for one pumping cycle. It is pro- 
portional to the dissipation ratio and reaches its maximum at 
times where the energy splitting reaches its minimum. 



was expectecP^, a zero temperature environment stabi- 
lizes the adiabatic theorem, i.e., does not influence the 
pumped charge. 

The solid red lines in Fig. [5] (a) represent the calcu- 
lation with perfectly closable SQUIDs, i.e., Ef^'' = 0, 
whereas the dotted black lines where calculated for real- 
istic SQUIDs with EY^'/EJ^'' = 0.03. Since both curves 
almost perfectly match, the pumped charge is hardly in- 



fluenced by the residual Josephson energy iJ™'" of the 
SQUID. Increasing the coupling to the environment sup- 
ports the charge transfer which is in agreement with the 
conclusions of Ref. HH Fig. [5] (b) depicts the total charge 
Qtot(''') transferred in one complete pumping cycle. As 
was expected, even for zero temperature environment the 
total charge is not quantized. It shows a similar behavior 
as the pumped charge when varying in temperature and 
coupling strength. Comparing both figures, the charge 
transferred due to super current only is also affected by 
varying temperature or coupling strength. The decrease 
of transferred charge for higher temperatures is due to a 
higher population of the excited state, which is responsi- 
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(b) 



FIG. 5. (a): Pumped charge Qp through the island as a 
function of couphng strength EqJo- The charge is plotted 
for different temperatures kBT/(2Tih) — 0,1,2,3 GHz (from 
top to bottom). The solid red lines represent calculations for 
perfectly closable SQUIDs, i.e. _E™'" = 0, whereas the dashed 
black lines were made for realistic SQUIDs. The weaker the 
coupling to the environment, the more the pumped charge 
tends to decrease. In Fig. (b) the total charge is depicted as 
a function of the coupling strength for different temperatures 
kBT/{2Tvh) = 0, 1, 2, 3 GHz (from top to bottom). 



ble for a back flowing current. Conclusively it should be 
mentioned that the obtained results are in a very good 
agreement with the results in Ref. [TW^. However, in [HI 
these results could only be obtained by including the non- 
secular terms, while our novel method using FCS allows 
for the use of the RWA. 



VI. CONCLUSIONS 

We explored an alternative way to calculate the 
amount of charge transferred in Cooper pair pumps. It 
turned out that introducing the counting fields allows 
us to perform the usual RWA and use the master equa- 
tion of the Lindblad form. We conjecture that FCS is 
in general the instrument of choice for charge calcula- 
tion whereas charge calculations via the current operator 
require a much more careful treatment. Our results are 
obtained numerically at this stage. A general analytic ar- 
gument for why FCS in combination with RWA produces 
accurate results, will be the subject of future work. 



Appendix A: Derivation of the Master equation 

As the system evolves periodic in time, we use the ba- 
sis of modified adiabatic Floquet statesSl |^„(i)). For 
a strictly adiabatic evolution, the system stays in one 
of its adiabatic eigenstates \n{to)). The modified Flo- 
quet states originate from the adiabatic eigenstates but 
include non-adiabatic corrections (for details se^^^ . 

Our starting point is a Markovian master equation in 
the interaction pictur^l^l 



|pPM(t) 



dsTvB[Hj{t), [Hj(t-s),ppM{t)(g)pB]] . 



(Al) 



To describe the time evolution of the system, we intro- 
duce the time evolution operator U — |^'„(t))((/)„(0)|, 
where is the time dependent state of the system. 

While using the adiabatic Floquet states, we rewrite the 
time evolution operator as 



U 



(A2) 



where En (t) are the instantaneous energies of the system 
and is the geometric phase. With the help of eq. 
( A2 ) we are able to perform a Fourier series expansion 



for {(l3n{t)\Ap\(j)„i{t)) as 



Ap{t) 



), (A3) 



with fl = 2tt/t. The Fourier coefficients are given by 



Ann',k = - [ dte-'^'"''{Mt)\Ap\<l>n'{t)). (A4) 



and the transition frequencies read 
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COnn'it) = E,At) - En{t) + {Ql - &h)/T. (A5) 



Eq. (Al), thus, becomes 
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dt 



PPM{t)= ^ 



nn' ^jj',k,l 

X (^\(l>n){'Pn'\PPM{t)\(l)f){(j)j\ - \(l)n){<Pn'\<Pj'){<Pj\pPM{t) 
X \(pj){(l)j>\ppM{t)\(l)n'){(pn\ - PPM{t)\(l>j){4>3'\4>n'){4>n 



(A6) 



with + Ofc) being the half-side Fourier (Laplace) 

transform of the bath correlation functions. These are 



given explicitly by 

poo 

Jo 

POO 

Jo 

(A7) 

The approximation made in the second step assumes that 
the eigenenergies of the system do not vary drastically on 
timescales of the bath correlation time. For large values 
of fifc sa ujnm the Fourier coefficients vanish, as the states 
\4>nit)) slowly oscillate. We arrive at 



n,j^k,l 

)\pT''-{t)) 

n,j,k,l 

-A:„^MA*j^MSjnP''p-''-{m 

X \c|>j{Vn^)){Mv^m)\Pp'"''''-{t)\M'p'J){M'PL)\-SjnA^^ 

(A8) 



Assuming adiabatic evolution, only small enough val- 
ues of flk contribute. Therefore the half-side Fourier 



transform of the bath correlation function can be ap- 
proximated as r{ujjni'Pm) + ^k) « r{ujjn{(Pm))- Using 
T{u>nm = 1/27 ('^nm) + we arrive at 
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Pp + iPp 



LS 



LS 



7(0) 



pI-^'-l','l'' 



^jkPp ^jk ^ i^jkPp ^jk 



(A9) 
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